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ABSTRACT 

We present a modification of the standard halo model with the goal of providing an im- 
proved description of galaxy clustering. Recent surveys, like the Sloan Digital Sky Sur- 
vey (SDSS) and the Anglo- Australian Two-degree survey (2dF), have shown that there 
seems to be a correlation between the clustering of galaxies and their properties such as 
metallicity and star formation rate, which are believed to be environment-dependent. 
This environmental dependence is not included in the standard halo model where the 
host halo mass is the only variable specifying galaxy properties. In our approach, the 
halo properties i.e., the concentration, and the Halo Occupation Distribution -HOD- 
prescription, will not only depend on the halo mass (like in the standard halo model) 
but also on the halo environment. We examine how different environmental depen- 
dence of halo concentration and HOD prescription affect the correlation function. We 
see that at the level of dark matter, the concentration of haloes affects moderately the 
dark matter correlation function only at small scales. However the galaxy correlation 
function is extremely sensitive to the HOD details, even when only the HOD of a small 
fraction of haloes is modified. 

Key words: cosmology: theory - cosmology: cosmological parameters - cosmology: 
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1 INTRODUCTION 

The modern language to describe analytically the clustering of galaxies is the halo model (e.g., [Cooray fc Sheth| | |2002[ ) for 
a review). In its original formulation, the halo model describes non-linear clustering of dark matter and can be applied to 
describe also the clustering properties of galaxies. 

In particular, the halo-model can be calibrated to describe galaxy clustering properties in several different ways, depending 
on how the observed properties of galaxies are to be related to the underlying dark matter halo - the so-called halo occupation 
distribution (HOD): using galaxy abundance, their spatial distribution via the two-point correlation function or the luminosity 
dependence as function of the halo mass. In the halo occupation distribution, it is customary to classify galaxies as central or 
satellite. These two classes of galaxies are then assigned different occupation numbers in the dark haloes. In this description, all 
observational properties are completely specified by the dark matter halo mass. The original halo model has been remarkably 
successful at describing the first moment statistics of the clustering of galaxies. 

In reality, however, galaxies are not easily divided in central or satellite, and the physical characteristics of galaxies of 
each type cannot be determined solely by the mass of the host halo: galaxies are more complex systems. Environment must 
play an important role in the process of galaxy formation, the most striking observational evidence being that clusters today 
have a much higher fraction of early-type galaxies than is found in the field. It has been known for more than three decades 
that there is a relation between galaxy morphology and density of the local environment starting from the results of 



Davis 



fc Geller| (|1976|); |Dressler| (|1980|l; |Postman fc Gellerj (|1984[), until most recent results e.g., |Hogg et alT] (|2002||; |Zehavi et a! 
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(20101. It is not clear if this effect can be completely ascribed to the fact that the most massive haloes are naturally found in 
overdense regions, and that most massive haloes form on average earlier, or if there is some extra environmental dependence; 
i.e., physical mechanisms such as ram-pressure stripping, harassment etc. ( Gunn fc Gott|19~7 2 Moor e~et al.|1996 19981 shape 
the properties of galaxies and operate in dense environments. We know that, although on average the star formation and 
mctallicity history are determined by the mass, their correlation properties arc not (Shcth et al. 2006 Mateus, Jimenez, & 
Gaztanaga 2008). In particular, the clustering of the properties of galaxies has been shown to depend on more parameters 



than mass: large-scale tidal fields have been shown to alter galactic spins ( Jimenez et al. |2010 1 and ellipticity ( Mandelbaum 
et al.|2006| > but can also alter galaxy clustering properties ( |Hirata|2009| [Krause fc Hirata|2010[ ). 



The fact that clustering properties of galaxies depend on galaxy internal properties, can be understood if dark matter 
haloes with different properties and formation histories, cluster differently and host different galaxy populations. Consider for 
example the so-called halo assembly bias; e.g. |Gao et al.| (j2005) found that the amplitude of the correlation function depends 
on halo formation time, thus haloes assembled at high redshift are more clustered than those formed recently. Models of 
galaxy clustering statistics make some simplifying assumptions; mainly that i) number and properties of galaxies populating 
a dark matter halo depend only on the mass of the host halo, and ii) clustering properties of dark matter haloes are only a 
function of their mass and not of the larger environment (this assumption is at the basis of the excursion set-formalism) . 

It is well known that the formation time of haloes depends on the halo mass and that objects formed at early times tend, 
on average, to be more concentrated than objects that formed recently. In the hierarchical galaxy formation model there is a 
correlation between galaxy type and environment induced by the fact that the mass function of dark haloes in dense regions 
in the Universe is predicted to be top-heavy. In a simple improvement of the classic halo model ( Abbas fc Sheth|2005| [2006 ) 
the most massive haloes only populate the densest regions and a correlation between halo abundance and environment is 
introduced. This approach matches the prediction of the excursion set approach; all correlation with environment comes from 



But there are indications that there could be a more complex dependence on the environment which is not fully described 
by halo mass e.g., |Sheth fc Tormen| ( [2004l ) and discussion above. It has been known for a while that the star formation history is 



perhaps the most affected quantity by environment, this is the so-called downsizing effect, e.g. (Cowie & Bargcr 2008 Heavens 



et al.|[2004 |. Furthermore, the recent zCOSMOS survey has shown how other variables are also affected by environment, in 
particular the shape of galaxy stellar mass function (Tasca et al.||2009| |Bolzonella et al.||2009|), and the formation of red 



galaxies at the low mass end ( Zucca et al.||2009 Cucciati et al.||2010 Galaz et al.||201o| . All these new observational results 



seem to indicate that the physical properties of galaxies are, at least in part determined by the environment in which they 
live, and that this dependence goes beyond the fact that the halo mass function is influenced by the local background density. 
Therefore, analyses of future surveys with a similar sensitivity as the zCOSMOS and, especially, large redshift coverage, may 
benefit from theoretical tools that can include environmental dependence. In practice the definition of "environment" may not 
be easy or univocal, and the definition may change as a function of the galaxy property under consideration. Here we simply 
set up the mathematical description assuming that "environment" has been already defined. We extend the standard halo 
model to include an extra dependence on the environment, classified as "cluster/node", "filament" and "void" although the 
treatment is general enough that other choices or interpretations of environment are possible. Motivated by the excursion-set 
model expectations, the extra dependence comes through the halo density profile and mass function for the dark matter and 
also through the HOD parameters for galaxies. In the excursion-set model this is given by the formation time but in our 
model we introduce an extra parameter set by the environment, leaving us one more degree of freedom to tune the HOD. As a 
starting point, here we introduce the formalism and show the effect of each of the parameters of the model on the correlation 
function. However we expect that the full potential of the model and the relevance of its features will appear when considering 



clustering statistics beyond the 2-point function of the over-density field such as e.g. the marked correlation functions (Skibba 



et al. 


2006 


(2011 


1. 



Sheth et al. 2006 Mateus, Jimenez, & Gaztanaga 20081. We defer this to forthcoming work Gil-Marin et al. 



This paper is organised as follow: in §2 we present the theoretical formulation of our extension of the model. We first 
introduce all the model parameters and ingredients, then in §2.4| we present the expression for the correlation function of dark 
matter and galaxies. In §3 we analyse how an scenario environments with different properties can affect to the correlation 
function. In §4 we present the conclusions of our work. Accordingly in the Appendices we review the basics of the standard 
halo model and the halo occupation distribution and give details of equations presented in §2. 



2 THE MODEL 

The halo model provides a physically-motivated way to estimate the two-point (and higher-order) correlation function of 
dark matter and galaxy density field. We review the classic halo model in Appendix A. Despite its simplicity, the model is 
extremely successful: the model's predictions have been compared with both simulations and observations, reproducing well 
the clustering properties. Nevertheless, it presents some limitations and shortcomings which we discuss next. The background 
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material presented in Appendix A is useful to set up the stage of motivating our extension of the model and to define symbols 
and nomenclature used. We thus refer the reader to Appendix A for definitions of many of the symbols used. 

In the halo model every halo is characterised only by its mass. This is correct at the leading order, and has been so far 
sufficient. Much more accurate measurements of dark matter and galaxies clustering will be available in the near future (e.g., 
SDSSIII, EUCLID, JDEM, DES, BigBOS^j] PannStarr, LSSlQ etc.) and a more sophisticated modelling may be needed 
to describe such data-sets. In fact there are indications already that the properties of galaxies depend somewhat on the 
environment, and not only on the host halo mass. See §1. 

It is clear that the mass has to be the primary variable in any halo model formalism. However in order to include 
some environmental dependence another variable needs to be introduced. At the level of the galaxy correlation function this 
"variable" can be modulating the HOD: in some environments galaxies can populate haloes in a different way than in other. 
At the level of dark matter we can introduce an environmental dependence through the concentration of the density profile. 
It has been observed that the relation between the mass and the concentration has a very high dispersion: the distribution 
around the mean concentration is approximately independent of halo mass, and is well approximated by a lognormal with rms 



o~inc = 0.3 (e.g., Sheth & Tormen ( 2004 1) . In the halo model formulation a relation between concentration c and mass m is 



assumed (see Eq. 33 1, however this high dispersion may indicate that the concentration depends on extra (hidden) variables, 



such as halo formation history (haloes which form at high redshift have a higher concentration), tidal forces etc., that in the 



end would depend on the environment. A possible approach to this problem was presented by Giocoli et al. (20101. It consists 
in treating the concentration as a stochastic variable and not only integrate over the mass but over all possible values of the 
concentration. In this work we will instead consider that the relation between concentration and mass may depend on the 
environment. 

In our model, we assume that the Universe contains three kind of structures or environments: nodes, filaments and voids. 
Halos lie either in a node or filament regions. Voids contain no haloes but occupy a large fraction of the Universe's volume, 
especially at z ~ 0. 

At the level of dark matter correlation function, we distinguish between haloes in the node regions (hereafter node- 
like haloes) and haloes in the filament regions (hereafter filament-like haloes) through the concentration of their profile. We 
assume that the node-like haloes have a profile with a concentration function c no d(m) and filament-like haloes with Cfu(m). For 
simplicity, in presenting our equations and our figures, we assume that these concentrations are constants, but the formalism 
can straightforwardly be generalised to allow these concentrations to be functions of the mass, as it is done in the halo model. 
The range of these values is expected to be between ~ 1 and ~ 15; these values (or functions) could be calibrated by N-body 
simulations. 

At the level of galaxy correlation function, we can decide to populate haloes in different ways, according to the type 
region where they lie. For instance, node-like haloes may have less satellite galaxies, or more likely to have a central galaxy 
than filament-like haloes. 

In this section we show how we can modify mathematically the halo model in order to take into account all these effects. 



2.1 Mass function 

As in the standard halo model (see Appendix A for details) , the halo number density is given by the mass function, which can 
be computed analytically in the extended Press Schechter approach or calibrated on N-body simulations. For this application 



we will adopt the Sheth & Tormen (1999) mass function. More massive haloes are more likely to be found in high-density 
environments, thus in principle the mass function can be made to be dependent on the environment and the dependence could 
be calibrated on N-body simulations. At this stage, this procedure is similar to the approach presented by [Abbas fc Sheth| 



(2005 20061. For this application the environment-dependence of the mass function is accounted for through the fact that 
the mass function depends on the mean matter density: p m , which we consider to be dependent on the local environment. 
However, this model is general enough to use any mass function with any environmental dependence. Let us consider a region 
i with a volume Vi and with a mean matter density p mi ■ If we consider this region of the Universe alone, then the number of 
objects of mass m in this i regions per unit of volume (of this region) and per unit of mass, namely m(m), is related to the 



mass function of the whole Universe as ( Abbas & Sheth 2005 ) 



n i (m) = [l + b(m)(y i -l)]n(m) (1) 

where ^ is the ratio between the mean matter densities of the region i and of the whole Universe: 3^ = Pmi/pm- It is also 
useful to define the volume fraction of this region i, Vi and the volume of the whole (observable) Universe V as Xi = Vi/V . 



See 

See http://www.lsst.org/lsst for more details 



http://bigboss.lbl.gov/index.html and http://arxiv.org/abs/0904.0468 for mote details 
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x y xy 



Nodes 3.8 x 10~ 3 73 0.2774 
Filaments 0.1368 5.28 0.7223 
Voids ~ 0.86 



Table 1. Volume fraction X and mean density fraction y. Values for node-, filament-like regions and voids used in this work. Data from 
( |Arag6n-Calvo et al.|2010^ 



Conservation of mass, imposes that the consistency relation must be satisfied: 

' x t y, = 1 



(2) 



where the summation runs over all kinds of structures, in our case nodes, filaments and voids. Since we assume that there are 
no haloes in voids (y v = 0), this equation reads: X„y n + Xfyj = l[^]ln other words, since the mass function is proportional 
to the mean matter density (see Eq. 19 1, if we consider a region with a mean matter density times denser than the mean 
density of the Universe, then the mass function of this region will be ^ times the one of the whole Universe (see Eq. [T]) 

For illustration, here the adopted values for these parameters are obtained from Aragon-Calvo et al. (20101 and are listed 
in table [J] In this work we choose these numbers as a fiducial values, but other values are possible. In fact, the values of 
these parameters will depend mainly on the definition of environment, and on the definition of what a node and a filament is. 

The filaments regions are more abundant than the node regions, but less dense. The voids occupy almost 90% of the 
volume. 



2.2 Halo density profile 

In the halo model, the halo density profile can be written as a function of mass and concentration: p(r\m, c). Usually the 
NFW profile ( Navarro et al.| 1996, 1997) is adopted, with an empirical relation between the concentration and the mass such 
as Eq. [33j For this application we will consider that these two variables are independent and that the concentration is set by 
the environment of the halo. Thus, we have two different profiles depending on whether the halo lies in a node-like region, 
p(r\m, c no d) or in a filament-like one, p(r\m, Cfu). In principle c„ oc j and c/y can be functions of the mass, which could be 
calibrated from N-body simulations. Here however, for simplicity, we will consider that these two variables are just constants. 



2.3 Halo Occupation Distribution 



In our approach, the standard HOD (described in Appendix A j ]6.6[ ) must also be modified to include a possible dependence 
on environment. 



As an example see the discussion in Zehavi et al. ( 2005 1: blue and red galaxies populate haloes of same mass in a different 



way (see Appendix A S 6. 6 and Fig. 16 1; blue galaxies tend to occupy low-density regions while red galaxies the high-density 
ones. 

Here we are interested in modelling generic properties of galaxies (star formation, metallicity etc.) not just colours. 
Whatever the property under consideration is, in our modelling galaxies can still be divided in two classes (nodes and 
filaments); for simplicity here we still call them "red" and "blue", but one should bear in mind that the argument is much 
more general. Thus we make the HOD to depend not only on the mass of the host halo but also on its environment (node or 
filament) by adopting two different prescriptions for nodes or filaments: the parameters M m i„, Mi and a of Eq. 42 and 43 



must be specified for node- and filament-like haloes. Let us continue with the working example of red and blue galaxies. While 
in the standard halo model, the red and blue galaxies are uniformly distributed inside haloes, in the extreme case of this 
environmental dependence, blue galaxies only populate haloes which are in filament regions while red galaxies only populate 
node-like haloes. However this is an extreme segregation and in a more realistic scenario there will be some mixing: node-like 
haloes host some blue galaxies and filament-like haloes can host some red ones. In order to describe this, we introduce the 
segregation index 5. When 5 = 0, there is no segregation, i.e. red and blue galaxies are distributed equally in filament and 
node haloes: this corresponds to the standard approach with no environmental dependence. In the other limit, if S — 1, 
blue galaxies live in filament haloes and red ones in node haloes: this is the extreme case of our extended halo model with 
environmental dependence. 



3 Here the sub-indices n, / and v accounts for the regions nodes, filaments and voids respectively 

4 In their work, the authors split their haloes in four types: clusters, filaments, walls and voids. Here we adopt the node values of X and 
y as their clusters haloes. For the voids we adopt X but we set y to 0, and we combine their filaments and walls values of X and y into 
what we call filaments. 
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Thus our adopted HOD equations are (see Appendix A for definition of HOD): 

9nod(m) = gred(m) + ^ ffMtte(m) 

g)a(m) = ^ 1 3 - \ 9biue ( m ) + 1 2 ~ gred (w) ■ (3) 

where 5 is the segregation index, g pr operty (m) denotes the average number of galaxies of a certain property in a halo of mass 
m and gj(m) denote the average number of galaxies in nodes (j — > nod) and in filaments (j — > /iZ). This <;(m) function is given 



by Eq. 42 and 43 for central (i — » cen) and satellite galaxies (i — > sat) respectively. These equations have 3 parameters: the 
minimum mass for a halo to host one (central) galaxy, M min ; the average mass for a halo to host its first satellite galaxy, M\; 
the power-law slope of the satellite mean occupation function, a. Taking into account that the HOD for red and blue galaxies 
may be different there are 6 HOD parameters. Moreover, if the two populations are not completely segregate but there is a 
degree of mixing, we have an extra parameter (see Appendix A for more details). In what follows, for some applications we 
will set g no d = gfU (i.e. g re d = guue) but in general the HOD parameters for g re d will be different from those of gu-ue and 
thus, in general, g nod / gfti . 

In next sections we will see how having two different HOD for nodes and filament regions modifies the correlation function 
respect to the case of having the same HOD for all galaxies. 



2.4 Two-point correlation function 

Following the above assumptions, we can compute the two-point correlation function for dark matter and galaxies. Since we 
have two types of structures (nodes and filaments), the correlation function is split in three terms instead of two as in the 
halo model|f] 

CM = e hlV (r) + f hlv (r) + e h2v (r) (4) 

where now £ lhlr, (r) accounts for particles in the same halo, ^ 2hlv (r) for particles in different haloes but of the same kind and 
£ 2h2r > (r) for particles in different haloes and of different kind separated by a distance r. In Fig. [I] we show schematically the 
contribution of these different terms. 



2-4-1 Dark Matter correlation function 

For dark matter particles the correlation function reads, 

On (r) = iZ " ( r ) + M + C h J V (0 (5) 
where each terms is, (see Appendix B for derivation) 

Cdm"(r) = V / dmX t m2n ^ m) f d 3 xu(x|m,cX|x + r||m, Cl ) (6) 

' ' Pm JVi 



i=l 
2 



~*dm 



t\ A j , v m'm{m') f „ m"m(m") f 3 , , ,. , f 3 „ . „. „ , 
(r) = y I am Xi ^ — - / am Xi ^ / ax u(x \m ,Cj) / a x u(x \m , CjJ 

- , J Pm J Pm Jy. J y. 



x £hh(W -x" +r\\m',m") (7) 

e 2h2 V , \ f , i v m'm(m') f „ m"n 2 (m") f 3 , , ,. , . f 3 „ . „, „ . 
Cdm ( r ) ~ 2 / dm X± ^ — - / dm X2 ^ / ax u(x \m ,ci) / ax u(x \m ,c 2 ) x 

J Pm J Pm Jy 1 Jy 2 

x Cwiflx' - x" + r\\m',m") , (8) 

where i = 1 refers to nodes and i — 2 to filaments. Here n(m) is the halo mass function, u(x\m,c) is the normalised density 
profile of the halo of mass m and concentration parameter c at a radial distance x (see Eq.|31[), p m is the mean matter density, 



£hh(d|mi, 7712) is the two-point correlation function of two haloes of masses mi and m 2 separated by a distance d (see Eq. 23 



for the model used here) and Xi and ^ are the volume and density fractions defined in E 2. 1 The integration J v d 3 x runs 
over all haloes' volume and j dm runs over all halo mass range (see Appendix A for computational details). In Fig[2]we show 
the contribution of the different terms. Here we have assumed c no d = 10 and Cfu = 2, just as an example to introduce the 
environmental dependence. According to the adopted values of X and y, the effect of the nodes is dominant at small scales 
(r < 0.3 Mpc/h) (red line in left panel). At intermediate scales (0.3 < r < 6 Mpc/h) the filaments dominate and at large 
scales (r > QMpc/h) both filaments and the cross term play an important role in the total correlation function (blue and 



5 The notation in this paper is the following: In means that both particles belong to the same kind of halo whereas 2n to different one. 
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Figure 1. In our modified halo model, dark matter haloes are of two kinds: node-like (red circles) and filament-like (blue circles). Because 
of this, our correlation function is split in three terms: the one-halo-one-?7 term (top-left picture) describes the interaction of particles 
inside the same halo, the two-halo-one-»7 (top-right picture) term involves particles which are in different haloes of the same kind. Finally 
the two-halo-two-r; term (bottom picture) considers the contribution of particles which are in different haloes of different kinds. 




Figure 2. Different contributions to the dark matter correlation function according to our model (Eqs.[6]-[8]|. Left panel: total correlation 
function Sl^m (black solid line), node contribution (red line), filament contribution rf* (blue line), cross contribution (green 

line) and linear dark matter contribution £; in (black-dotted line). Right panel: total correlation function £,dm (black solid line), 1-halo 
term S}^^ (red line), 2-halo-same-halo term f j 1 *' (blue line), 2-halo-difierent-halo term £^ 2?) (green line) and £,n n (black-dotted line). 



green lines in left panel). This is expected due to the fact that the filaments are more abundant than nodes but nodes are 
more concentrate than filaments. 

In the right panel we show the contribution of the different terms of Eq. [5] As we expected, the lftljj-term (red line) term 
dominates at sub-halo scales whereas 2hl-q-term (blue line) dominates at large scales. The cross correlation term between 
nodes and filaments, 2/i2?7-term (green line) is less important than filaments at large scales but is not negligible. In the left 
panel, the shape of the filament and the node contribution is set by the value of the concentration we have chosen. We 
will come back to this point in the next section to see how the concentration parameter affects the shape of the correlation 
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function. In both panels the linear correlation function is also plotted (black-dotted line). The ratio between £u n and at 
large scales can be defined as an effective large scale dark matter bias, 

f , mn(m) , , . 

b dm = / dm } > b(m) = 1 (9) 
J Ph 



by definition of b(m) (see Eq. 25 1 this dark matter bias is one if we integrate over all range of masses. Also the effective large 



scale bias can be defined for the node and filament contribution, this effective bias is given by 

b dm = Xi [ dm mniijn) b(m) = XM . (10) 

J Pm 

Here i stands for either nodes or filaments and the relation bdm = b™^ + 6^ = 1 is satisfied. Note that the only difference 
between the the dark matter bias of nodes and filaments is due to Xi and 3^- We will see that if we rescale the mean density 
p m to the mean density of nodes or filaments (Xip~ mi ), then the effective bias for nodes and filaments is the same. 

Note that when we refer to the node and filament contribution to the total correlation function, we refer to the terms of 
the sum of Eq. [5] These terms are different from the correlation function we would obtain if we only took into account the 
node- or filament-like haloes (we call it 'pure' node and filament contribution). In the first case, each term is divided by the 
total dark halo-matter density squared, p^, while in the second ('pure') case instead it would be divided by the density of 
node- or filament-like haloes^] In Fig.|4](left panel) the contribution of these 'pure' terms is shown for the dark matter. Note 
that the shape is the same as in Fig. Inland the only difference is that the normalisation is shifted by a factor (A'Jj) -2 . 

2.4-2 Galaxy correlation function 

In this case the environmental dependence can be introduced not only through the concentration parameter of the host halo, 
but also through the HOD: galaxies do not populate in the same way filament- and node- like haloes even if the host halo 
has the same mass. In other words, if we have two galaxy populations (red and blue) with different HOD, e.g., red galaxies 
are more abundant in node haloes and blue in filament haloes, then an environmental dependence arises. This is the case 
presented here. 

As before, the galaxy correlation function reads, 

^(r) = tiller) + e g h J V (r) + C( r ) (") 

where each term is, 

tlal V ( T ) = ^2 [ dmX t n !^ ri ^ 2gi en (m)g s l at (m)u{r\m,c l ) + (g- at {m)) 2 [ d 3 yiu{x\m, d) u(\yL + r\m,d) (12) 
i=1 J n gai L J J 

/-2hln / \ \ ^ f j / , // v 2 n i{ m )ll/i(m ) r cen , /. C en . //. . , / //, cen / tt\ sat/ /\ 

igai ( r ) = 2.^ \ dm dm Xi — 2 [g t {m)g l (m )£ hh {r\m ,m ) + 2 g t (m )pi (m ) x 

i=i J n gai 

x J d 3 x u(x'\m\ci)£hh(\it + r\\m' ,m'') + gi at (m')g'i at (m") J d 3 x' 'd 3 x" u(x'\m ,Ci)u(x"\m" ,Ci) (13) 
x £«i(l x - x + rj|m , m J 

*2h2-n i \ r- f j J j // v v n\{m )n2(m ) r cen , cen/ tl\s- / i / //\ . cen/ l\ sat/ ff\ 

C 9 ai ( r ) = 2 dm dm X X X 2 — 2 [g ± (m)g 2 (m )£ hh (r\m ,m ) + g 1 (m)g 2 (m ) x 

J n gal 

x J d A x" u(x"\m" ,c 2 )^hh(W + r\\m ,m") + g2 en {m")g'l at (m) J d 3 x' u(x'\m, ci)^ hh (|x' + r| \m , m") (14) 

. sat/ l\ sat/ ll\ I i3 / j3 // / 1 1 I \ / II II \ /- /\ I II . II / //\ 

+ 9i ( m )92 \ m ) / d x d x u{x\m,c\)u[x \m , c 2 )t i j l i l (\x — x +r\\m,m ) 
Again i — 1 refers to nodes and i = 2 to filaments; g J (m) is the average number of galaxies that lie in a halo of mass m, the 



superindex j denote the type of galaxy: central (j — > cen) or satellite (j — > sat) (see S6.6 for details); n ga i is the mean number 
density of galaxies, i.e. the total number of galaxies divided by the total volume, 

2 

n ga i = ^2Xin aa u (15) 
i=i 

and n ga i j is the mean number density of galaxies inside haloes of type i, i.e. the total number of galaxies inside haloes of type 
i divided by the volume these haloes occupy, 

n ga ii= J dm ni(m) [gi en (m) + g s t at (m)] . (16) 

6 In order to obtain these 'pure' node and filament terms we have to multiply the mean density p m by the factor Xyi, in order to 
obtain the mean density of nodes (i = 1) or filaments (i = 2) 
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Figure 3. Same notation as in Fig. j2j Left panel: £ sa ; (black solid line), £™°f (red line), (blue line) and Cg^f' 7 (green line). Right 

(green line). In both panels 5dm (black dashed line) and £iin 



plot: £ ga [ (black solid line), £ 9a ; ^ (red line), C gai ^ (blue line) and 
(black dotted line). 



c 2h2-q 



gal 



As before we can define an effective large scale galaxy bias as, 

S^Xif dm 7 ^ (gr(m) + gr\m)) b(m) . (17) 



'Jgal 



In this case the effective large scale galaxy bias for nodes and filaments is given by, 



Kal = %i I dm n% }" b ' {gT n {m) + gt at (m)) b(m) . (18) 

J flgal 

As before b ga i = b^f + b„^ v Note that for the galaxy bias, the difference between the node and the filament bias not only 



depends on the terms Xi and as in the dark matter bias (Eq. 10 1, but also on the HOD of each population. This means 
that even rescaling the mean number of galaxies n ga i to the mean number of galaxies of nodes or filaments, the biases will be 
different as we will see in the next section. The behaviour of these terms is illustr ated in Fig. [3] For simplicity ci = C2 = c(m) 
given by Eq. 33 is adopted, along with the HOD for galaxies introduced by Zehavi et al. (20051: log 10 M m in = 12.72, 



log 10 Mi = 14.08 (the masses are in M Q /h), a = 1.37, / cen = 0.71, f$ at = 0.88, a c A T = 0.30 and < = 1.70 (see Appendix A 
f 6.6 for detailed definition of these HOD parameters) with maximum segregation: 5 = 1. This means that blue galaxies only 
populate filaments and red galaxies only populate nodes. The left panel shows the correlation functions of nodes (red line), 
filaments (blue line), cross term (green line) and the total (black solid line). In the right panel, the different lines correspond 
to the terms of Eq. |11| \h\r\ (red line), 2h\r\ (blue line), 2K2r\ (green line) and the total contribution (black solid line). For 
comparison, the dark matter correlation function (black-dashed line) and the linear power spectrum (black-dotted line) are 
also shown. 

In the left panel of Fig. [3] we show the effect on the total correlation function (black solid line) of red galaxies (red line), 
blue galaxies (blue line) and the cross term (green line). We see that according to this HOD, red galaxies dominate the 1-halo 
term (r < 4 Mpc/h), whereas all three terms contribute to the 2h term, being the cross term the most important. 

We can also compute the 'pure' node and filament galaxy terms, as we did in the dark matter case. These terms are 
shown in Fig. [4] (central panel). As before, the shape is the same in both Figures but in the second one, the lines normalisation 
is offset by (<Vj34) -2 - I n the case of the 2h2r\ term, the rescaling goes as ^ ai n — > £g a fV( 2<-fi3^i X2 3^2) in order to obtain cross 
bias at large scales defined as b cross = ^/b nod bf l1 . This is due to the mean density definition as it is explained below. 

In Fig. [4] we show the 'pure' contribution of the node and filament terms for dark matter (left panel) and galaxies (central 
panel). The correlation function for nodes/red galaxies is considerably higher than that for filaments/blue galaxies at all scales. 
This is just an effect of rescaling the mean density/number of galaxies from the total mass/galaxies over the total volume 
to the mass/galaxies in filaments or nodes over the total volume. This rescaling goes like for the mean density /number 
of galaxies and like (A 1 ;^)" 2 for £. Since at small scales £ ~ Ai3^i, the ratio between the node and the filament 'pure' £ 
terms is just (Xfyf)/(X n y n ) ~ 2.5 and the pure correlation function for nodes is 2.5 times higher than that of filaments 
at small scales in the case of dark matter and a bit different in the case of galaxies due to the different HOD, as it can be 
seen in Fig. [4] (left and central panel). However at large scales £ ~ (Xiyi) 2 , which compensates for the rescaling of the mean 
density /number of galaxies. Thus, in the case of dark matter both filament and node pure correlation function are exactly 
the same since both effective large scale bias are the same after the rescaling (see Eq. 10 1. However in the case of the galaxy 
correlation function, the effective large scale bias is considerably different for nodes and filaments due to they have different 



© 0000 RAS, MNRAS 000, 000-000 



Halo Model and Environment 9 




Figure 4. Separate contribution for node- (red lines), filament- like (blue lines) haloes, the cross term (green line) for dark matter (left 
panel) and galaxies (central panel). The ratio between galaxies and dark matter correlation function is shown in right panel. Black lines 
correspond to the total contribution (node+filaments+cross). 



HOD (see Eq. |18[ ), being higher for red galaxies than for blue ones. In the right panel the ratio £ ga i/£,dm is shown for nodes 
(red line), for filaments (blue line), for the cross term (green line) and for the whole sample (black line). Both at small and 
large scales the ratio is higher for node-like haloes. 



3 ANALYSIS 

In this section we perform a more exhaustive analysis to how HOD parameters, segregation and concentration affect to the 
dark matter and galaxy correlation function. We analyse how every parameter of the HOD affect the dark matter and galaxy 
correlation function while the other parameters are kept fixed. We also analyse how the concentration parameter and the 
segregation affect to the both dark matter and galaxy correlation function while any other parameter is fixed. To gain insight 
we first start with the standard halo model (i.e. 5 = and only one HOD prescription for al haloes) and then we move to our 
extended halo model with environmental dependence. 



3.1 Halo Model 



1.2 



Here we adopt Eqs 
and a) affect £ ga i- Here M, 



and 43 as a parametrisation of the HOD and we analyse how the different parameters (M m in, Mi 
sets the minimum mass for a halo to have galaxies, Mi is the mass of a halo that on average 

for 



hosts one satellite galaxy and a is the power-law slope of the satellite mean occupation function (see Appendix A S 6.6 



more details). We set the parameters to fiducial values which are very close to the ones proposed by Kravtsov et al. (20041: 
M m i n — 1O 11 M0, Mi = 22 x 10 11 Mq and a = 1.00. In the following section we analyse how changing each one of these 
parameters affects the galaxy correlation function and the ratio £gai/£dm- Since filament galaxies are the most abundant, 
the effects we find in this case would be qualitatively very similar to the effect of changing the HOD of filament population, 
keeping the nodes HOD fixed at the fiducial model. 



3.1.1 Central Galaxies 

We start by varying M m i„, the minimum mass of a halo to host at least one galaxy, keeping all the other HOD parameters 
fixed at their fiducial values. In Fig.[5]the different solid lines correspond to different values of M min : 10 11 (red line), 5 x 10 11 
(blue line), 10 12 (green line), 5 x 10 12 (pink line) and 1O 13 M0 (orange line). The left panel is the galaxy correlation function 
and the right panel is £ ga l/£dm- For reference, in the left panel, the dark matter correlation function £d m is also plotted 
(black-dashed line). 

We can see that at large scales the galaxy correlation function is just a boost of the dark matter correlation function, 
whereas at small scales the shape is modified. The enhancement of the amplitude of the correlation function at large scales can 
be explained just as a redistribution of galaxies: removing galaxies from low mass haloes (M < M m i n ) causes a large-scale bias 



increase because the function 6(m) (Eq. 25 1 increases with the mass. Therefore at these large scales we expect that only the 
amplitude of the galaxy correlation function changes with M m in (and not the shape that is given by £«„). On the other hand 
at small scales the 1-halo term of the galaxy correlation function becomes important; this term does not depend explicitly 
on £dm, therefore the shape of the correlation function does not have to be the same. The behaviour of Fig. [5] is qualitatively 



very similar to the findings of Berlind & Weinberg ( 2002 ) 
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Figure 5. Effect of changing M m i n in the HOD model. In the right panel the galaxy correlation function £ ga i is shown, whereas in the 
left panel the ratio (, ga i/£,dm ls shown. The different colours correspond to different M m i n values: 10 11 Mq (red line), 5 X 10 11 Mq (blue 
line), 10 12 M Q (green line), 5 X 10 12 Mq (pink line) and 10 13 Mq ( orange line). In the left panel, black-dashed line corresponds to the 
dark matter correlation function £dm- 




Distance (Mpc/h) Distance (Mpc/h) 



Figure 6. Effects of changing Mi in the halo model. As in Fig.[5]the left panel is £ ga i and the right panel £, ga l/£,dm- Different colours 
are different values of Mi: 22 X 10 10 Mq (red line), 44 X 10 10 Mq (purple line), 22 X 10 11 Mq (blue line), 22 X 10 12 Mq (green line) 
and 22 X 10 13 Mq (orange line). In left plot £<j m is a l so represented (black dashed line). 

3.1.2 Satellite Galaxies 

Here we keep fix M m in and a to the fiducial values 10 11 Mq and 1.00 respectively and we vary Mi from 22 X 10 10 Mq (red 
line) to 10 13 Afg (orange line). Recall that this is the minimum mass for a halo to host at least one satellite galaxy. 

In Fig. [5] we show the effect of varying this parameter. The left panel shows the galaxy correlation function and the right 
panel £ ga i/£dm for different values of Mi (see caption for details). In the left panel, ^ m is also plotted (black-dashed line). 
At sub-halo scales we observe that high values of Mi suppress the correlation function. This effect is expected: the higher 
Mi, the less satellites galaxies per halo and therefore the 1-halo term is suppressed. We also observe that for very high values 
of Mi (> 22 x 1O 12 M ) at very small scales £ ga i shows an inflection point around 0.3Mpc//i. This is due to the interaction 
between the central galaxies of different haloes. Since for these values of Mi there is (almost) no satellite galaxy per halo, 
we can observe the interaction between central galaxies of the smallest haloes (those of a mass ~ M m i„), which are the only 
ones that at that distances can contribute to the 2-halo term on those scales. Such feature is similar to the one of the 2-halo 
term of Fig. [3J just with the difference that in that case M m in ~ 10 12 ' 72 Mq, whereas now Mmin = 10 11 Mq and therefore 
such feature happens at shorter distances. 

At large scales we also see that the correlation decreases as we increase Mi. This is also expected: according to Eq.[l7|, low 
mass haloes are weighted more in contributing to the large scale bias, as one decreases Mi the large scale correlation function 
also decreases. In other words, for high values of Mi all haloes count the same (1 central galaxy per halo) independently of 
their mass. 
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Figure 7. Effects of changing a in the halo model. As in Fig. [5] left panel is £ ga ; and the right panel £ ga i/£dm- Different colours are 
different values of a: 0.80 (red line), 1.00 (blue line), 1.10 (green line), 1.30 (pink line), 1.40 (orange line) and 1.50 (purple line). In the 
left plot £,dm is the black-dashed plot. 

Now we consider the variation of a from 0.8 to 1.4, keeping M m in — 10 11 M@ and Mi = 22 x 10 11 M©. This parameter 
indicates how the number of satellite galaxies increase with the mass of the halo. A steeper slope a indicates more extra 
satellites per halo mass increment. 

In Fig. [7] we show the effect of changing a in the galaxy correlation function (left panel) and in £, ga i/£,dm (right panel). 
The different colours indicate different values for a (see caption for details) and the black-dashed line is the dark matter 
correlation function. 

As expected, when a increases, the galaxy correlation function also increases, especially at sub-halo scales. At large scales 
the effect of increasing a is the same as increasing Mmin'- it boosts the dark matter correlation function without changing its 
shape. This is because increasing a we weight more the more massive haloes and therefore the bias increases. At small scales 
there is also an enhancement of the correlation function for the same reason, but in this case the shape of (, ga i is qualitatively 
different from the shape of £dm, because as in the case of changing M min , there is not direct relation between the 1-halo terms 
of £, ga i and £dm- This behaviour was also noted by |Berlind fc Weinberg] ( [2002] ) . 

3.2 Extended Halo Model 

In this section we use again the HOD parametrised by Eq. [3] using as g cen (m) and g sat (m) for filament- and node- like haloes 
Eqs. [42] and [43] but with different parameters. In particular we keep fixed the HOD parameters for filament haloes to the 
fiducial values: M m %n = 10 n M©, Ml = 22 x 1O 11 M0, a — 1 and we allow to vary M m in, Mi and a only for red galaxies. We 
assume a concentration c(m) given by Eq. [33] and begin with a maximum segregation index 5 = 1 (this will later be relaxed). 
Since the filaments are the most common structure, changing their HOD would give a similar effect as changing the HOD of 
the whole population, that we studied in the previous section. Therefore, in this section we focus on how changing the HOD 
of the less abundant population (in this case the nodes) can affect to the total correlation function. 

3.2.1 Central Galaxies in nodes 

We set Mi = 22 x 1O 11 M and a = 1.0 for node-like satellites and we vary M min from 10 11 M Q to 10 13 M Q . 

In Fig. [8] we show the effect of changing M m in for node-like galaxies. The different colour lines are £ ga i for different values 
of Mmin for node-like galaxies (see caption for details) and the black-dashed line is £dm- In t ne left panel we show the total 
correlation function. We see that the total galaxy correlation function is not very sensitive to this parameter. We observe only 
a small to moderate change, much less that that observed when we changed M m in of the total population of haloes. This 
can be due to the fact that if we increase M m in less node-like haloes are populated, so there are less red galaxies and the 
total correlation function is dominated by the filament-like galaxies. In the right plot we show the corresponding effect on the 

£,gal / ^dm • 

3.2.2 Satellite Galaxies in nodes 

Here we keep M m in = W 11 Mq and a — 1.0 fixed and we allow Mi to vary. In Fig. [9] we show the effect of changing Mi in 
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Figure 8. Effect of changing M m ; n for node-like haloes in the extended halo model. Left panel is £ ga ;; right panel is £ ga i/£,dm- The 
different colours arc different values for M min : 10 11 M (red line), 5 X 10 11 M (blue line), 10 12 Mq (green line), 5 X 10 12 Mq (pink 
line) and 10 13 Mq (orange line). In the left panel §<j„i is also plotted (black-dashed line). 




Figure 9. Effect of changing M\ in node-like haloes. As in Fig.|8]in left panel £ ga i is shown and in the right panel £, ga l/£.dm is plotted. 
The different colours correspond to different values for Mi: 22 X 10 10 Mq (red line), 44 X 10 10 Mq (purple line), 22 X 10 11 M (blue 
line), 22 X 10 12 M (green line) and 22 X 10 13 M (orange line). In the left panel, the black-dashed line represents £,dm- 



the galaxy correlation function (left panel) and on ^ ga i/^dm (right panel). The different colour lines represent different values 
of Mi . As before, the black-dashed line is the dark matter correlation function. Recall that the effect of increasing Mi is to 
reduce the number of satellite galaxies for low-mass haloes in the node-like regions. We see that £ ga ; is only sensitive to the 
change of Mi in the range from 22 x 1O 1O M to 22 X 1O 11 M . For values larger than 22 x lO n M the correlation function 
saturates. In particular we see that £ ga i is especially sensitive to Mi at small scales where the 1-halo term dominates. This 
means that the number of satellites in the low-mass and node-like haloes plays an important role in the total correlation 
function. 



In Fig. 10 we show the effect of changing a in node-like galaxies keeping fixed M min and A/i to 1O 11 M and 22 x 1O 11 M 
respectively. We explore the regime from a = 0.80 (red line) to a = 1.50 (purple line) for the £ 9a ; (left panel) and for (, ga i/(,dm 
(right panel). As before, £dm is the black-dashed line. We observe the same effect observed in Fig. [7] As in the Mi case, £ S<1 ; 
is especially sensitive to a at small distances where the 1-halo term dominates. This may indicate that for these values of 
X and y the satellite galaxies in the node-like haloes have an important role at small scales on the total galaxy correlation 
function in spite of not being the dominant population. 



3.3 Concentration 



In Fig. 11 we show the effect of changing the concentration parameter for both dark matter (left panel) and galaxy correlation 
function (central panel) . Also in the right panel we show the effect on £ ga i / (,dm ■ The different lines are for different concentration 
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Figure 10. Effect of changing a on nodes-like haloes in our halo model. As in Fig. [8] in left panel £ ga ; is shown whereas in the right 
panel £,g a i/£,dm is plotted. The different colours show different values of a: 0.80 (red line), 1.00 (blue line), 1.10 (green line), 1.30 (pink 
line), 1.40 (orange line) and 1.50 (purple line). In the left plot, the black-dashed line represents §dm- 
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Figure 11. Correlation function for different values of the concentration parameters: c no( j = 10 & c*n = 2 (red line), c no( j = 2 & 
Cfn = 10 (blue line), c not j = 10 & Cfu = 10 (green line) and c not j = 2 & Cfu = 2 (orange line). The black-dashed line represents the case 
for c(m) according to Eq. |33| In the left panel, the dark matter correlation function is plotted, in the central panel the galaxy correlation 
function and in the right plot 5ga//5dm- 

values of nodes and filaments: c no d = 10 & Cfu = 2 (red line), c no d = 2 & Cfu = 10 (blue line), c no d — 10 & Cfu = 10 (green 
line) and c no£ j = 2 & Cfu =2 (orange line). For comparison we also have plotted c(m) according to Eq. 33 (black dashed 
line). Here we have fixed the HOD to the one of Zehavi et al. (20051: log 10 M m in = 12.72, log 10 Mi = 14.08 (the masses are in 



Mq /h) , a = 1.37, /o en = 0.71, fo = 0.88, a^" = 0.30 and crjf = 1.70 (see Appendix A ije^for detailed definition of these 
HOD parameters) and adopted the maximum value for the segregation index, 5 = 1. Note that significant changes appear 
at small scales (r < 1 Mpc/h). Also these changes are only important for the dark matter correlation function. On the other 
hand, for galaxies are almost negligible. 

In general the effect of changing the concentration can be summarised up as a change in the slope of the the small scale 
correlation function (especially the dark matter one), with a more concentrated profile (green line) corresponding to a steeper 
correlation function. 



3.4 Segregation 

In Fig. [12] we explore the effect of having the two galaxy populations mixed or segregated. For two population of galaxies 
completely segregated the segregation index is 5 = 1 (green line); when the two population are perfectly mixed this segregation 
index is S = (red line) (see Eq. [3] for definition of 5). Intermediate cases 1>5>0(5 = 0.5 for blue line) describe a certain 
degree of mixing between the two populations, but always with red galaxies being the most abundant population in node-like 
haloes and blue galaxies more abundant one in filament-like haloes. The galaxy correlation function and £ iga i/£,dm are plotted 
in the left and right panel respectively. In the left panel, the black-dashed lines is the dark matter correlation function, ^ m 
We observe that at large scales, the more mixed the two population of galaxies are, the higher is the correlation function. This 
can be understood from the fact the red galaxies are more biased (see Fig. |4|. If the two population are completely mixed 
(5 = 0) all haloes will contain red galaxies and therefore, in total, there will be more red galaxies than if only the node-like 
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Figure 12. Effect of the segregation parameter. Left panel is £ ga i and right panel £ ga i/tidm- For 5 = 1 (green line) the red and blue 
galaxies are completely separated: red galaxies live in node-like haloes and blue galaxies in filament-like haloes. For 5 = (red line) 
the red and blue population are perfectly mixed in node- and filament-like haloes. For S = 0.5 (blue line) an intermediate situation 
is presented: red galaxies are more abundant in node like haloes and blue galaxies in filament-like haloes, but some degree of mixing 
between these two population is allowed. In the left panel, the black-dashed line represents §<j m . 



haloes could host this red population (as it happens with 5 = 1). Therefore with increasing 5 the relative importance of the 
more biased galaxy population decreases, thus decreasing the total correlation function. However this effect it is very small 
in the total correlation function. 



4 DISCUSSION & CONCLUSIONS 

In the classic halo model the host halo mass is the only variable specifying halo and galaxy properties. This model has been 
remarkably successful at describing the first moment statistics of the clustering of galaxies. However, environment must play 
an important role in the process of galaxy formation, the most striking observational evidence being that clusters today have 
a much higher fraction of early type galaxies than is found in the field. In addition, when looking at statistics beyond the 
two-point density correlation function, there are indications that the simplest halo-model may be incomplete. 

When looking at the two-point galaxy correlation function for current, low-redshift surveys, the data can be well modelled 



without such an additional dependence on local density (beyond that introduced by the host halo mass, see e.g., (Skibba & 
Shcth 2009)). However, the statistical power and redshift coverage of forthcoming surveys implies that an extended modelling 
may be needed. 

In this work we have presented a natural extension of the halo model that allows us to introduce an environmental 
dependence based on whether a halo lives in a high-mass density region, namely a node-like region, or in a low-mass region, 
namely a filament-like region. At the level of dark matter, the secondary variable (the environment-dependent variable beyond 
the shape of the mass function) is the concentration of the density profile of the haloes: haloes which live in nodes regions 
may have a different concentration than haloes which live in filaments regions. According to this idea we present the dark 
matter correlation function, £dm, of this new extended halo model: Eqs.[6j[8] In the classic halo model the correlation function 
is usually split in two terms: the one-halo and the two-halo term. In our extended halo model we have 3 terms: the one-halo 
term and 2 two-halo terms, depending whether the two particles belong to haloes of the same or different type (see Fig. [I] for 
a graphical example). We have explored the contribution of these terms to the dark matter correlation function and we have 
found that the contribution of particles that belong to haloes of the same and different kind is similar (see right panel of Fig. 
[2J. We have also seen that according to the values of table [TJ the contribution of the node-like haloes starts to be dominant 
at small scales and is subdominant at large scales (see left panel of Fig. [2j). We also have analysed how different values of the 
concentration affect the correlation function. At the level of dark matter, changing the concentration parameter between 2 



and 10 only affects moderately the correlation function at small scales (see left panel of Fig. 11 1. 

We also have extended our environment-dependent halo model to galaxies, computing the galaxy correlation function, 
(, ga i: Eqs. |12fl4| This has been done by modulating the Halo Occupation Distribution recipe with the environment. In the 
galaxy correlation function, the HOD is our secondary variable rather than halo concentration. This means that according to 
our model, node-like haloes host galaxies in a different way than filament-like galaxies. In this paper we use a simple model for 
the HOD of only 3 variables: the minimum mass of a halo to host its first (central) galaxy, M m i„; the mass of the halo to have 
on average its first satellite galaxy, Mi and the way how satellite galaxies increase with the halo mass, a. Thus we have two 
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different HOD, one for red-galaxies and another for blue-galaxies. Here we have chosen the colour as an example, although 
another property besides the colour, like star formation rate or metallicity may be used instead. Of course more sophisticated 
HOD could also be used, but starting with a simple prescription helps gaining physical insight. We analyse how changing the 
different parameters of the HOD affects to the galaxy correlation function (see Fig. 5)10 1. We see that, even when we change 
the parameters of the HOD of the minority population (in this case the node-like haloes), there is a considerable change in 
the £ g ai, especially for the parameters Mi and a, i.e., for the satellite population. Therefore changing the HOD of only a 
small fraction of haloes affects considerably to the total galaxy correlation function. We also have explored how the mixing 
or segregation between the two population of red and blue galaxies may affect to £ ga i- We expect that the size of this effect 
will depend on the HOD of each one of these populations: the more different are the HOD of red and blue galaxies, the more 
the effect of segregation. 

It is reasonable to expect the dependence of galaxy properties and clustering on environment to be very complex in 
details. But, being a second-order effect, it is likely that a simplified description that captures the main trends will be all that 
in needed in practice. The model presented here is simplified in several ways: 1) the environment is only divided into nodes, 
filaments, and voids and it does not have a continuous distribution; 2) that dark halo clustering only depends on mass not on 
environment. The formalism introduced here allows one to implement 2) straightforwardly. 

A natural extension of this model is to make a continuous dependence on the environment, instead of splitting it in nodes, 
filaments and voids. However this extension presents several issues. First of all we would need a continuous dependence on the 
environment of the mass function and also of the HOD, which is not yet available. Secondly, nowadays N-body simulations 
start to provide information of haloes according a discrete number of environments: this discretization could be taken as a 
first approximation to the continuous distribution. While this avenue is worth pursuing, introducing a more complex model 
based on a continuous environment dependence goes beyond the scope of this paper. 

We envision that the extension of the halo model presented here will be useful for future analysis of large scale structure 
of the Universe, especially analyses that account for the physical properties of galaxies. Current photometric and spectroscopic 
large-scale surveys are beginning to gather not only the position of galaxies but also some of their physical properties, like 
colour, star formation, age or metallicity, and with future surveys the accuracy of these measurements will increase. We expect 
that the physical properties of galaxies depend on the environment, making these datasets the suitable ground to apply the 
extended halo model. 

The halo model presented here uses analytic expressions for the mass function, halo density profile, and also it uses a 
given HOD for a certain magnitude-selected galaxies. We understand that these analytic functions may be inappropriate for 
comparison with data. In the present paper we wanted to have an expression for the dark matter non-linear power spectrum, 
that is fast and easy to compute and a reasonable approximation to reality; we however were interested in the relative effect 
of including environment dependence and less concerned with absolute accuracy of the fit. In a forthcoming paper we are 
planning to apply this technique to fit observational data, in this case we would have to further refine our prescription to 
compute ^dm (see Smith et al. (20031 for a numerical approach to the halo model description). This is still work in progress 
but we believe it goes beyond the scope of the present paper. 

The environmental halo model presented here is especially suited to model the marked (or weighed) correlation formalism 
| |Harker et al.||2006] [Skibba et al.|2006) , consisting on weighting each galaxy according to a physical property and removing 
the dependence of clustering on local number-counts. Here we have laid the foundations for modelling a survey's marked 
correlation, which treatment will be presented in a forthcoming work. 
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6 APPENDIX A. THE HALO MODEL 



In this section we review the basics of the standard halo model. This background material set up the stage for motivating and 
introducing our extension of the standard halo model and will define symbols and nomenclature used. 



The halo model pioneered in Peacock & Smith (20001; Ma & Fry (20001; Seljak (20001 then thoroughly reviewed in 



Cooray & Sheth (20021, assumes that all the mass in the Universe is embedded into units, which are called dark matter 



haloes or simply haloes. These haloes are small compared to the typical distance between them (non-linear evolution makes 
the evolved Universe dominated by voids). For this reason, the clustering properties of the mass density field, <5(x), on small 
scales are determined by the spatial distribution inside the dark matter haloes, and the way they are organised in the space 
is not important. On the other hand the statistics of the large-scale distribution is not affected by the matter distribution 
inside haloes but only by their spatial distribution. In this work we assume all the time halo exclusion, treating haloes as hard 
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Figure 13. Left panel: global mass function (solid line), mass function in node-environment (dashed line) and mass function in filament 
environments (dotted line). Right panel: Ratio of the mass functions of the left panel. n no d(m) /n(m) (dashed line), nf i i(m)/n(m) 
(dotted line) and n no d(m)/nfu(m) (solid line). 



spheres (sharp cutoff at the virial radius). In the case that two particles that belong to different haloes, are separated by less 
distance that the sum of the virial radii of their host halos, we avoid counting them in the computation of the correlation 
function. 



6.1 Mass function 



The number density of collapsed haloes of mass m per unit of mass at a given redshift z, n(m,z), can be computed using 
the Press-Schechter formalism ( Press fc Schechter|1974 1 identifying the present collapsed haloes with the peaks of an initially 
Gaussian field, 



n(m, z) 



2p-n 



da(m) 



dm 



(19) 



ma(m) 

where p m is the mean density of matter in the Univers^J o(m) is the rms of the power spectrum linearly extrapolated 
to z — filtered with a top-hat sphere of mass m and v = 8^ c (z)/a 2 (m). Here, S sc (z) is the linearly extrapolated critical 
density required for spherical collapse at z, and is given by 8 3C (z) — 1.686/_D(z), where D(z) is the growth factor. Following 
this formalism, structures with a linearly evolved density fluctuation higher than this threshold value, will collapse. The ST 
approach (Sh eth fc Tormen|1999 1, based on ellipsoidal collapse, predicts a mass function of the form, 

f{v)u = A{p) (1 + (,„)-*) (g) 1/2 exp {-qu/2) (20) 

with p = 0.3 and q = 0.75 and the normalisation factor A(p) = (1 + 2~ p r(l/2 — p)/r(l/2)) _1 . Note that by construction the 
matter density is given by, 



p™. = 



dm n(m, 



(21) 



In order to introduce an environmental dependence in the mass function, one could think of rescaling Eq. [19] to account the 



differences of densities. However Mo & White (1996) noted that dark halo abundance in dense and underdense regions do not 



differ by just a factor like this. A more complex model was proposed by | Abbas fc Shet h (2005), 

n(m;Mi,Vi) = [l + b(m)6i]n(m) (22) 

where, n(m; Mi, Vi) is the mass function of a region of volume Vi which contains a mass Mi, b(m) is the bias of a halo of mass 
m (see Eq. [25] for details), and Si is defined as Mi/Vi = p m (l + Si) In Fig. 13 (left panel) we have plotted the global mass 
function (Eq. |19|) (solid line), a high-density environment mass function (dashed line) and a medium-density environment 



mass function (dotted line) using in all cases the ST mass function (Eq. 201. For these two last mass functions we have used 
respectively the relative density parameters of table [T] corresponding to nodes and filaments. On the right panel we show the 
ratios of these three different mass functions: n no d(m)/n(m) (dashed line), nfu(m)/n(m) (dotted line) and n no d(m) / n fu(m) 



7 According to a ACDM flat universe (Q m = 0.27, = 0.73 and h = 0.7) the matter density is given by p m = — where 
Ho = 100?t. Using the fiducial cosmology assumed here this value IS pm — 7.438 x lO lo M //i(Mpc/?t)" 3 
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(solid line). We see that the ratio of the number of massive to low mass haloes is larger in dense regions (nodes) than in less 
dense regions (filaments) since the slope of the solid curve of the right panel of Fig. 13 is positive. 



6.2 Bias 

For large separations the bias is the relation between the correlation function of two dark-matter haloes of masses m! and 
m" separated by a distance r at a given redshift z, £,hh{r, z, m! , m"), and the underlying dark matter linear power spectrum, 
£lin(r, z). While computing the exact form of £hh(r, z, m' , m") is a somewhat delicate matter, excellent results are obtained 
by using 

£hh(r, z,m , m") ~ b(m' , z)b(m" , z)£i in (r, z) (23) 

where b(m,z) denotes the large-scale linear bias factor for haloes of mass m at redshift z. While using as a proxy for 
£hh is strictly incorrect because there may be mildly non-linear contributions and because at small separations haloes are 
spatially exclusive, on these scales the signal is dominated by the one-halo term. While the fit to N-body simulations can be 
further improved by using e.g., higher-order perturbation-theory prediction instead of e.g., Smith et al. (20081, we will 
not pursue this here as it is beyond the scope of the present paper. 

To be consistent we should use a bias prescription derived from the extended Press-Schechter formalism using the peak 
background split. To the lowest order the bias is, 

91n[n(m, z')] 



b(m, z) — 1 



dSsc(z') 



The bias of an object m at time z is given, at lowest order of 5, by (Cole & Kaiser 1989; M o fc White|199"6 1 



b(m, z) 



1 + 



D[z) 



Ssc(z) _ 1 

cr 2 (m) S BC (z) 



(24) 



(25) 



where, D(z) is the linear growth factor, S sc (z) is the critical threshold for collapse and is given by S sc ~ 1.686/D(z), a(m) is 
the rms of the power spectrum linearly extrapolated at z = filtered with a top-hat sphere of mass m and q is the parameter 
introduced in the Sheth 



Tormen mass function (see Eq. 20 1. This formula has been confirmed by N-body simulations giving 

i 

, we 



an excellent agreement. Although the complete formula should have the extra term, 2p/(5 ac (z)D(z)) ^ 
have checked that the effect of this term is about 1% and can be safely neglected for our application. 



1 _i_ ( Aw 

1 + { CT (m)2 



6.3 Density profile 

According to our definition, a halo is a set of particles that forms a gravitationally bound and thermodynamically stable 
system. Therefore, we consider that this system satisfies the virial theorem, i.e. a halo is virialised by definition. From the 



spherical collapse model (Gunn & Gott 19721, a halo is considered to be formed and virialised when its density reaches a 

187r 2 p„ 



certain threshold value: p V i r . Here we will adopt p v i 
spherical collapse model. The corresponding virial radius is then, 



which is the typical value that can be found using the 



1/3 



(26) 



m 

2n \3p m 

In the halo model, this will be the size of an halo of mass m. 

On the other hand, cosmological simulations have shown that the density profile inside isolated haloes follows a universal 



profile given by (Navarro et al. 111996 19971, 



P(r\r s ,p s ) 



Ps 



(27) 



{r/r a )(l+r/r s ) 2 

where r s is the scale radius of the halo and p s /4 is the halo density at scale radius. However it is more natural to work with 
the concentration c, and the mass of the halo m, instead of p s and r s when we describe the halo profile. These variables are 
defined as, 



; Tvir I T s 



and the mass of the halo has to satisfy 



From this last equation, we can write, 



r r v, 
Jo 



r(m) 



drAwr p(r\r a , p s 



(28) 



(29) 



Ps = 



47T7-3 



ln(l + c) - 



1 + c 



(30) 
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Figure 14. In the halo model, the correlation function is split in two terms. The one-halo term (left) describes the clustering of particles 
inside the same halo. On the other hand, the two-halo term (right) involves only particles of different haloes. 



We prefer to use normalised profile of a halo u as a function of distance from its centre for given halo virial mass and 
concentration parameter, defined by, 

p(r\m,c) 



which satisfies the condition 



u(r\m, c) 



1 = 



d x u(x\m, c) . 



(31) 
(32) 



6.4 Concentration 

In principle the halo density profile depends on two independent parameters: either p s and r s , or c and m. However N-body 



simulations indicate that there is a relation between the concentration and the mass. Seljak (20001 propose 

n -0.2 

W III 

c(m, z) — 



9 



m 



(33) 



where m* is the mass of a typical collapse halo at z = 0. There are many other parametrisations of the relation between the 



concentration and the mass (e.g., Neto et al. (20071; Zhao et al. ( 2009 1) . However we do not expect that the results of this 
work to depend significantly on this choice. Unless otherwise stated, for the present application we thus adopt Eq. [33] for the 
relation between concentration and mass. 

This relation goes in the direction one may have expected: less massive haloes on average form earlier, when the Universe 
is more concentrated and therefore have a higher concentration than more massive ones. One should bear in mind however 
that there is a large dispersion around this mean relation, which may correspond to a "hidden parameter" such as local 
environmental effects, tidal effects, mergers etc. 



6.5 Two-point correlation function 

In the halo model, the two-point correlation function of dark matter particles contained in haloes is defined as, 

Cd m (r) =< S(x)<5(x + r) > (34) 

where <5(x) = p(x)/p m — 1. This function can be split into 2 terms, depending on whether the two particles at distance r 
belong or not to the same halo: 

Cdm (r) = f dm (r ) + Cdm (r ) • (35) 

Cdm( r ) i s called the one-halo term and accounts for particles in the same halo; dm( r ) i s called the two-halo term and accounts 
for particles belonging to different haloes. Therefore the properties of the mass density on small scales are described by £d m (r"), 
whereas on large scales are given by £d m (r*). A graphical description is shown in Fig. 14 
According to the halo model formalism, these terms are, 

,h , . f , m 2 n(m, z) f l3 , \ , nr .\ 

t,dm( r )Z) — / dm — _ 2 — - / d xu(x|m) u(|x + r\\m) (36) 

J Pm\Z) Jv 

& 2h f \ f j ,m'n(m',z) f „m"n(m",z) f 3 , , f 3 „ „ „ , „ , „ 

£ dm (r, z) — dm — - — p-r — / dm — / d x ulx \m ) / ax u(x \m Khh\ x x + r , z, m , m ) (37) 

J Pm{Z) J Pm{Z) J v J v 
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1000 




0,1 1 10 
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Figure 15. Total dark matter correlation function £<jm (solid line), one-halo term (dashed line), two-halo term £3™ (dot-dashed 
line) and linear power spectrum £; in (dotted line). 



The integration J v d 3 x means over all haloes' volume, and J dm runs over all mass range. For computational reasons we have 
to adopt some limits in the mass integral. Setting m ma x = 10 Mq is enough as far as the value of the integral do not change 
if we increase this limit. If we set the minimum mass limit to 10 9 Mq it is also enough for the 1-halo term, because the less 
massive haloes contribute at very small distances (less than 0.1 Mpc/h). However the contribution of low mass haloes to the 
2-halo term its not negligible. In order to compute these integrals we use the method described in Yoo et al. ( 2006 1 which 
consists in breaking the 2-halo integral in two parts and approximate the low mass halos as points without inner structure, 



1 mn(m,z) , 
dm f -b(m, z) 

Pm 



d x.u(x\m) 



dm 



mn(m, z) 

pm 



b(m, z) / d 3 xu(x\m) + 1—1 
Jv . Jm 



mn(m, z) 

pm 



b(m, 



(38) 



where we have used the approximation of Eq. [23] to explicit the mass dependence of ^hh- In particular we have set mth to 
lO 9 A/0 and we have checked that reducing this value do not affect the result of the integral 
In Fig. 



15 



we can see the contribution of these two terms. £<im( r ) dominates at scales typically smaller than the virial 
radius, whereas £j^(r) does it for larger scales following the shape of £; in slightly shifted by the effect of the bias (see Eq. 23 1. 



6.6 Halo occupation distribution 

The halo model provides a framework to also model galaxy clustering: the complicated galaxy formation physics would 
determine how many galaxies form in a halo and their sampling of the dark matter halo profile. Thus, the shape of the 
one-halo term would be modified accordingly. 

Typically the HOD assume a centre-satellite distribution of galaxies inside each halo. This means that a galaxy is placed 
at the centre of the halo and may be surrounded by satellite galaxies distributed according to some statistics. The average 
number of galaxies that lie in a halo of mass m is (e.g., |Sheth| ( [2005] ), 

g W {m) = J2Np(N\m) (39) 

N 

where p(N\m) is the probability density of N galaxies are formed in a halo of mass m. 

We can also define the nth factorial moment of the distribution p(N\m) of galaxies in haloes of mass m, 

<?(n)(m) = ]T N(N - 1) • • • (N - n + l) P (n\m). (40) 

JV 

Here we assume that p(n\m) follows a Poisson distribution, although other choices are also possible. Under this assumption, 
we can state that for satellite galaxies, 

9(n)(m) = (ff(i)(m)r ■ ( 41 ) 
In particular we will only be interested in the first moment gn)( m ) (the mean number of galaxies in a halo of mass m) and 
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in the second moment (j( 2 )(m) that we will treat as the square power of the former. Therefore hereafter we write g( 1 )(m) as 
g(m) to simplify the notation. 

The typical way to parametrise the mean number of centre galaxies is to think of the mean number of the central galaxies 
as a step function, 



C<:: 7) / \ 

g (m) = 



1 if m>M„ 
if m < M m 



The satellite galaxy distribution can be parametrised as a power law, with the same mass cut as central galaxies, 

sat [ \ I 7TI \ cen / \ 

9 ( m > = \T^r) 9 M 



(42) 



(43) 



where M m ; n sets the minimum mass for a halo to have galaxies; Mi is the mass of a halo that on average hosts one satellite 
galaxy; and a is the power-law slope of the satellite mean occupation function. All these parameters can be tuned to fit 
observations and in general depend on the type of galaxy under consideration. 



In the literature there are several ways to populate haloes with galaxies depending on their type: LRGs ( Reid & Spergel 



2009), field galaxies (Kravtsov et al. 20041 and many others (e.g., Zehavi et al. (2005)). On one hand, Kravtsov et al 



(2004) using high-resolution dissipationless simulations find that the parameters of the HOD according to Eqs. 
M min ~ 10 11 M e /h, Mi ~ 22M min and a ~ 1. On the other hand, 



Zehavi et al. 



42 



and 



43 



(2005) based on SDSS observation state that 



for a volume-limited sample with magnitude in r-band (M r ) < —21 the values of these parameters are 

Mi — io 14 08 Mq/H and a = 1.37 (see Fig. 17 left panel), although these values depend strongly on the magnitude selection 
criteria. Therefore it is clear that the parameters of the HOD depend strongly on the way galaxies are selected: whether one 
deals with a volume-limited sample of whether a magnitude selection criteria has been applied. It is interesting to note that 
some authors (e.g., Zehavi et al. ( 2005 1) introduce a colour dependence through the HOD: haloes of the same mass host 



different kind of galaxies, red and blue. Since blue galaxies tend to live in low-density regions (filamentary regions) and red 
galaxies tend to live in high-density regions (node regions), this can be though as different HOD parametrization for different 



kind of environments. In particular, Zehavi et al. (20051 introduce the fraction of blue galaxies in a halo of mass m, fb(m), 



as a function of the halo mass and also as a function of galaxy type: central or satellite. Since red galaxies are more common 
in high mass haloes the authors adopt a function for fb which is decreasing with halo mass. For central galaxies a log- normal 
function is adopted: 

(log M m 

in , 



b (m) = fo 

whereas for satellite a log-exponential function is used 

/sat [ \ rS 
b (m) = f 



cxp 



n cen'2 
M 



exp 



lo. 



•Sin m - 



log 10 M min 



T sat 



(44) 



(45) 



where m is the mass of the halo, M m in is the minimum mass of a halo to host one galaxy, fo is the value of the fraction 
for a halo with mass M m i n and au is a parameter that characterise how fast the blue fraction drops. Also these parameters 
depend strongly on the magnitude-selection criteria. For M r < —21 they are: log 10 M min — 12.72, f§ e " = 0.71, fo at = 0.88, 
a™ = 0.30 and er^"' = 1.70. The average number of blue and red galaxies becomes then, 



gliue{m) = f b {mYg l (m) 
gled{m) = [l - fbijn)^ g'(rn) 

where the average number of galaxies of type i, g l {m) is given by Eqs. 



42 



and 



43 



(46) 

according the values of the parameters 



used in Zehavi et al. (20051. In Fig. 16 we illustrate the corresponding HOD for Eq. 46 for red (red lines) and blue (blue 



lines) galaxies, and for the total number of galaxies (black line). In the left panel, we show the separate central and satellite 
contributions: the solid lines are central galaxies and dashed lines satellite galaxies. In the right panel we show the combined 
effect of the contributions. As stated by the authors, for haloes just above M m i n , blue central galaxies are more common. 
However above 2M m i n , central galaxies are predominantly red. Note that for blue galaxies there is a minimum in the total 
number of galaxies that occurs when a halo is too massive to have a blue central galaxy but not massive enough to host a 
significant population of blue satellite galaxies. As a general trend, we can say that according to this kind of HOD model: i) 
the lowest mass haloes have blue central galaxies, it) higher mass haloes have red central galaxies but with a significant blue 
fraction in their satellite population and finally in) the highest mass haloes host red central and satellite galaxies. |Zehavi "et] 



al. ( 2005 1 considered that the two populations are perfectly mixed and proceeded to compare the HOD predictions for the 



projected correlation function with SDSS data. This kind of difference between the HOD of blue and red galaxies motivates 
us to further explore the environmental dependence of the HODs of haloes of different density regions and in particular to 
consider that different environment can have different ratios of red and blue galaxies. 
Eqs. |36| |37|can then be modified to describe the clustering of galaxies: 
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Figure 16. HOD with colour dependence according to |Zehavi et al,| ( |2005| l. Left plot: Central and satellite galaxy contribution (solid and 
dashed lines respectively) for red and blue galaxies (red and blue lines respectively). Right panel: Total galaxy number for all galaxies 
(black line), red galaxies (red line) and blue galaxies (blue line). 




Log M h 



Figure 17. Left panel: HOD for the -21 > M r sample. All galaxies (solid line), central galaxy (dashed line) and satellite galaxies (dotted 
line). Centre panel: Dark matter correlation function §dm (solid line), galaxy correlation function £ ga ; (dashed line) according the HOD 
of the left panel and linear correlation function ^i™ (dotted line). Right panel: ratio between £ ga ; and £dm- 



Cgaiir,z) = i dm 



n(m, z) 
n 2 Az) 

qal V / 



2 g cen (m)g sat (m)u(r\rn) + (p sat (m)) 2 / d 3 x.u(x\m) u(|x + r\rn) 



(47) 



dm'dm" n{m '> Z Ji m "> z) 

n g al( Z ) 



cen / I \ cen / a \ j- / I u\ , r\ cen / i \ sat / f/\ i id / / 1 1 /'\ ,, 

g {rn)g (m )£ hh (r, z,m ,m ) + 2 g (m )g (m ) / d x u{x \m ) X 



s- n I i I ' ' f \ i sat/ f\ sat/ ff\ I j3 / t3 // / l\ l\ / ll\ ll\/. /\ I It . \ I ll\ 

x £hh(|x + r| , z,m , m ) + g (m )g (m ) / axax u(x \m )u(x \m )^hh(|x — x + r|, z, m , m ) 



/-2h / \ 



where the mean number of galaxies per halo is given by 

n ga i{z) = / dmn(m,z) (g ce "(m) + g sat (m)) 



(48) 



(49) 



In Fig. [17] (left panel) the |Zehavi et al.| ( |2005| ) HOD is plotted (M < -21 sample, log 10 M mi „ = 12.72, log 10 Mi = 14.08 
and a = 1.37). The dashed line corresponds to the contribution of central galaxies, the dotted line to satellite galaxies and the 
solid line to the total. In the central panel we show the correlation function for dark matter and galaxies: £ ga i (solid line) , ^ m 
(dashed line) and £; in (dotted line), according the same HOD. In the right panel, the ratio between £ ga i and is plotted. 
We see that at scales smaller than ~ 1 Mpc/h, when typically the one-halo term dominates, the ratio (, g ai/£dm is highly scale 
dependent, whereas at larger scales, then the two-halo terms dominates the ratio £, ga i/£,dm is scale independent. In particular 
it appears that at scales ~ 0.1 Mpc//i there is a minimum for t; ga i/idm with a value of 8. At large scales (> 10 Mpc/ft) £ ga i/£dm 
takes a value around 1.3. This is expected from the halo model formalism, where at small scales the 1-halo term dominates 
whereas at large scales the 2-halo term does it. From Eq. |47|48| it can be seen that at large scales the 2-halo term can be 
expressed as a function of whereas this does not happen for the 1-halo term. This yields us to the point that at large 
scales the ratio between galaxies and dark matter it is just a constant, whereas at small scales the relation it is much more 
complicated. 
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7 APPENDIX B 



Here we set the basics of the derivation of Eqs. |6|8| of our extended halo model. This derivation is very similar to the one of 
classic halo model presented for the first time by |Scherrer fc Be rtschinger ( 1991 1. 

We define the overdensity <5(x) as 5(x) = p(x)/p — 1. Then the dark matter correlation function is, 



£ dm (r) = <5(x)5(x + r)> = ^(p(x)p(x + p)> - 1 

P 



(50) 



Where (• • • } denotes the ensemble average. On the other hand, the density field at x is the sum of densities of node- and 
filament-like regions: p(x) 
halo of this class: 



p nod (x) +p /I( (x). We also can write the contribution of the nodes/filaments as the sum every 



/5 j (x) = S~] miu(yi — X||mi, Cj) (51) 

i 

where the summation takes place over all haloes of type j (nodes or filaments), m,i and Xi are the mass and the position (of 



the centre) of the ith halo and u is the normalised profile defined in Eq. 31 From Eqs. 50 and 51 it is clear that we can 



split £dm in 3 terms: Cdm* 1 anc ^ £ 2 ' 1 depending on whether the particles belong to the same halo and depending on the 
particles belong to the same kind of halo (in the case of the two-halo term). The derivation of each terms to Eqs.[6][8]is similar 
and are based on the introduction of the Dirac delta, 



— > / dm d 3 x 5% (x — Xj)<5p (m — m,) 



and on the definition of the mass functions in the node- and filament-like regions as, 

1 



n.,-(m) = 



Vi 



V"d 3 x5f,(x - Xi)5 (m - 



(52) 



(53) 



where the index j denotes the type of environment, either node or filament, and Vj the volume of this environment. 

Using these formulae and following the formalism presented by |Scherrer fc Be rtschinger] ( |199ij ), Eqs. |6)8| can be obtained 
after some algebra. 
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